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Abstract 

For ensemble of bosons trapped in a ID harmonic potential well we have found 
an analytical formula for the canonical partition function and shown that, 
for 100 trapped atoms, the discrepancy between the grand canonical and the 
canonical predictions for the condensate fraction reaches 10% in the vicinity of 
the Bose-Einstein threshold. This discrepancy decreases only logarithmically 
as the number of atoms increases. Furthermore we investigate numerically 
the case of a 3D "cigar-shape" trap in the range of parameters corresponding 
to current BEC experiments. 
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Recently, Bose-Einstein condensation (BEC) in trapped atomic gases |],|2| has been re- 
alized. The trapped atomic cloud possesses two remarkable features: First, the system is 
small enough so that finite particle effects are potentially observable, and second, particle 
interactions are weak. The thermodynamics of such a system is an interesting and rich area 
for scientific analysis. 

The equivalence of the grand canonical and fixed- iV canonical descriptions of a statistical 
system is an old question widely discussed in the textbooks on statistical mechanics H|J. 
For a bosonic gas, where grand canonical fluctuations of the ground state population become 
large at and below the Bose-Einstein threshold, such an equivalence is not obvious. It is 
shown that in the thermodynamic limit iV — > oo both ensembles give the same predictions 
for the mean values of occupation numbers even in the absence of particle interactions 
||. Furthermore it is well-known that for large N, interactions between particles lead to 
suppression of fluctuations in the grand canonical ensemble [|J . However, for a finite system 
with a mesoscopic number of particles, the equivalence of the two ensembles is not ensured. 

The main scaling laws for the fluctuations in an ideal canonical bose gas are derived by 
Fujiwara et al. [[J. The fixed- N bose statistics is shown to be closely related to Gentile's 
grand canonical intermediate statistics ||. Krauth has performed fixed- N finite temperature 
Monte-Carlo calculations for a 3D harmonic potential 0. Although the main subject of 
the paper is the role of interactions, it is shown also that for macroscopic numbers of 
particles, the noninteracting grand canonical and canonical ensembles agree very well. These 
conclusions are consistent with the numerical results of Politzer ||. In the present paper we 
consider mesoscopic values of number of particles (N ~ 100) confined in a one-dimensional 
harmonic trap || and in a three-dimensional "cigar-shape" trap. We show that the grand 
canonical/canonical deviations in this case are substantial. 

Consider an ensemble of N noninteracting bosons confined in a ID harmonic potential in 
thermal (but not in diffusive) equilibrium with a large reservoir. The population distribution 
among the different energy levels of the iV-particle system will be given by the Boltzmann 
law: 
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p([n])~exp{-/3£([n])} (1) 
[n] = {n ,ni, ...,n s , ...| ^n s = , (2) 

where = J2 s n s e s is the A^-particle energy for the given configuration of occupation 

numbers [n], e s = fckjs (s = 0, 1,2, ...) is the single particle energy spectrum, u is the har- 
monic oscillator frequency, j3 = 1/ksT, and T is the temperature of the system. Note that 
in the harmonic oscillator case, the iV-particle energy is quantized as 

E([n]) = %wK, (3) 

where K — J2 S n s s - 

To calculate mean occupation numbers of the oscillator states, we need to know the 
partition function Q and its derivatives. We show below that in the ID harmonic oscillator 
case, occupation numbers may be calculated analytically as finite sums of finite products. 
The canonical partition function Q(j3,N) can be represented by a power series of x = 
exp(—j3huj): 

oo 

Q(j3,N)= "£x K r(K,N). (4) 

K=0 

The micro canonical partition function 

T(K,N)= 1 

In] 

Sn s =N,Hn s s=K 
N 

= E E 1 (5) 

N'=0 [„]' 

T,'n s =N',T,'n s s=K 

equals the number of representations (partitions) of K as an unordered sum of at most iV pos- 
itive integers. Here N' is the total population of the excited states, [n] = {n 1; n2, n s , ...} 
is a particular configuration of excited state occupation numbers, and the "primed" sum 
£' = denotes a sum over the excited states. 

According to a well-known number theory theorem [JLO] the number T(K, N) of partitions 
of K with at most N parts equals the number V(K, N) of partitions of K with parts not 
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exceeding N. Hence, the canonical partition function (||) is nothing else but the generating 
function for the restricted partition function V(K, N) ||10||: 



Q(P,N) = J2x K nK,N) 

K=0 
N 1 



= n ■ (6) 

N=l L X 

Derivatives of the partition function Q s = —f3^ 1 {dQ/de s ) can not be found directly from 
the expression (§) which is specific for the ID harmonic oscillator. Instead, we have found 
a general recursion relation between the canonical partition function and its derivatives: 

Q s {(3, N + l)= exp(-(3e s ) (Q s (f3, N) + Q((3, N)) . (7) 

This relation can be applied to any fixed- N, noninteracting, bosonic system. 
Finally, the mean occupation numbers are given by 

N _ N 

N=l Ar=jv 

This expression is easy to analyze in the continuous limit with respect to N. For example, 
below the BEC threshold, the condensate population is approximately given by 

^ w i-^(io g (WM + c + o(i)) 

^1-^, (9) 

■L C 

where C ~ 0.5772 is the Euler constant. The transition temperature is given by 

k B T c ( const k B T c \ 

iV = ^7 lo H^^J' (10) 

where the choice of const is a matter of convention. Note that the thermodynamic limit (||) 
coincides with the one predicted for grand canonical statistics 0. 

Now we are ready to compare the canonical and grand canonical predictions for the con- 
densate population (n ). In Fig. 1 we plot the population of the ground state for different 



numbers of particles. For the grand canonical predictions we simply repeat the finite-system 
calculations of ||. To facilitate the comparison, we made the same choice const = 2 in 
expression ([T0|). Both curves approach the thermodynamic limit ([]) as the number of par- 
ticles increases. However, for a finite number of particles the discrepancy between the two 
models is quite significant. In the vicinity of the BEC threshold, the relative deviation 
((nQ T ' canon ') — (nQ anon ')) / {riQ anon ') decreases slowly with N and goes from 10% for 100 atoms 
to 5% for 10, 000 atoms. We have checked that this deviation decreases according to a 
1/ ln(iV) scaling law for a fixed T/T c . Note that the rate at which both the grand canonical 
|| and canonical @ populations approach the thermodynamic limit also obeys this law. 

We turn now to the 3D trap. To our knowledge there is no simple analytic expression for 
the canonical partition function in this case. Numerically, it can be calculated by integration 
of the grand canonical partition function in the complex plain of chemical potential |)]||. 
Indeed 

Q(P,N)= £ exp{-/?£([n])} 

[n] 

£n s =7V 

= J2 6 Xn a ,N exp{-(3E([n])} 

[n] 

= JL dfiexp(-Nti)Z(J3,n), (11) 

Zm J-ni 

where 

OO 1 

m»)= II i n vT ( 12 ) 

is the grand canonical partition function, uj a (a = x, y, z) are the trap frequencies, and the 
expression S qjq > = (27ri)~ 1 f^™ <i^exp[(g — q')£\ for the the Kronecker delta has been used. 
Derivatives of the partition function can be expressed through Z(/3,fi) in the same way. 

In Fig. 2 we plot the condensate fraction as a function of temperature for both grand 
canonical and canonical ensembles. We have chosen the "cigar-shape" configuration uj± = 
17.78 u> z , where u± = u x = u y . The three dimensional Bose-Einstein transition temperature 
is given by 



# = 0.(1)^ 



(k B T c 



fa Yia=x,y,z 

+ ^^W^ + 0(]fc B T c /M , (13) 

1 fa \.\a=x,y,z UJ a 

where the second line is the finite-iV correction ||. Here ga(z) = Y^?=\Z? 1 3 d is the Bose- 
Einstein function. For comparison, we have also plotted the thermodynamic limit 

( "°>-l-f£V. (14) 



N \T C 

For 100 particles, depending on the temperature, the system exhibits both 3D and ID 
characteristics. At T ~ 0.4 T c the temperature reaches the zero-point energy hu±/2 for 
tranverse oscillations. The grand canonical/canonical discrepancy is less than in the purely 
ID system but is still close to 10%. 

In the above discussion we neglected particle interactions. To estimate the importance 
of interactions in our model, we consider the "worst" case of zero temperature where the 
spatial density is the highest and therefore the interactions are strongest. For typical Ioffe- 
Pritchard trap parameters = 2tt x 101 Hz, u z = 2ir x 5.7 Hz, N = 100) for sodium 

atoms (scattering length a = 92 Bohr, atomic mass M = 23 amu) the mean-field corrections 
to the oscillation frequencies are quite small: 5ux_ = 0.03 u± and 5uj z = 0.15 u±. To estimate 
the corrections, we minimized the Gross-Pitaevskii energy functional with a ground state 
oscillator wave function seeded with unknown frequencies |12|] . Note that for the parameters 
chosen, the system exhibits a BEC transition at T c = 6.5 nK. 

We acknowledge fruitful discussions with H.D. Politzer, T.H. Bergeman, J.H. Thywissen, 
E. Heller, L. You, M. Prentiss and W. Ketterle. M.O. was supported by the National Science 
Foundation grant for light force dynamics ^PHY-93- 12572. C.H. was supported by Harvard 
University. This work was partially supported by the NSF through a grant for the Institute 
for Theoretical Atomic and Molecular Physics at Harvard University and the Smithsonian 
Astrophysical Observatory. 
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FIGURES 

FIG. 1. The condensate fraction for the ID harmonic oscillator as a function of temperature. 
Both grand canonical and canonical predictions are shown. The straight line is the thermodynamic 
limit (9). 

FIG. 2. The condensate fraction for a 3D "cigar-shape" trap as a function of 
temperature. Both grand canonical and canonical predictions are shown. Here 
uj ± = 2tt x 101 Hz, u z = 2ir x 5.7 Hz, N = 100, T c = Q.bnK. The thermodynamic limit N = oo 
(14) is also shown. 
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